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Coupled Variations of Fundamental Couplings and Primordial Nucleosynthesis 
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The effect of variations of the fundamental nuclear parameters on big-bang nucleosynthesis are 
modeled and discussed in detail taking into account the interrelations between the fundamental 
parameters arising in unified theories. Considering only 4 He, strong constraints on the variation of 
the neutron lifetime, neutron- proton mass difference are set. These constraints are then translated 
into constraints on the time variation of the Yukawa couplings and the fine structure constant. 
Furthermore, we show that a variation of the deuterium binding energy is able to reconcile the 
7 Li abundance deduced from the WMAP analysis with its spectroscopically determined value while 
maintaining concordance with D and 4 He. The analysis is strongly based on the potential model 
of Flambaum and Shuryak that relates the binding energy of the deuteron with the nucleon and a 
and ui mesons masses, however, we show that an alternative approach that consists of a pion mass 
dependence necessarily leads to equivalent conclusions. 
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I. INTRODUCTION 
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Primordial nucleosynthesis (BBN) is one of the most 
successful predictions of the standard hot big-bang 
model. Its success rests on the concordance between the 
observational determinations of the light element abun- 
dances of D and 4 He specifically, and their theoretically 
predicted abundances [3, 0| • Furthermore, measurements 
of the CMB anisotropies by WMAP [3[ have led to preci- 
sion determinations of the baryon density or equivalently 
the baryon-to-photon ratio, rj. As ij is the sole parameter 
of the standard model of BBN, it is possible to make very 
accurate predictions [1, H, d, 0] and hence infer the ex- 
pected theoretical abundances of all of the light elements 
(including 7 Li, 6 Li, 9 Be, 10 Be and n B). 

At present, a discrepancy between the predicted abun- 
dance of 7 Li and its spectroscopically determined abun- 
dance persists. The 7 Li abundance based on the WMAP 
baryon density is predicted to be [5j]: 



dances in these stars yields Q Li/H = (1.23^ 
and more recently Li/H = (1.26 ±0.26) x 10~ 10 . This 
value is in clear contradiction with most estimates of the 
primordial Li abundance, as also shown by Q who find 



7 Li/H 



4.26t°;g x 10 



-10 



(2) 



7 Li/H 



4 - 15 ±8:« x 10 



-10 



(1) 



The systems best suited for Li observations are metal- 
poor halo stars in our Galaxy. Analyses of the abun- 
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In both cases, the 7 Li abundance is a factor of ~ 3 higher 
than the value observed in most halo stars. 

There have been several attempts to account for the 
discrepancy between the BBN/ WMAP predicted value 
of 7 Li/H and its observational determination. These in- 
clude depletion mechanisms due to rotationally-induced 
mixing and/or diffusion. Current estimates for possible 
depletion factors are in the range ~ 0.2-0.4 dex [l(|. 
However, the negligible intrinsic spread in Li leads 
to the conclusion that depletion in halo stars is as low 
as 0.1 dex. It is also possible that the stellar parameters 
used to determine the Li abundance from the spectro- 
scopic measurements may be systematically off. Most 
important among these is the effective temperature as- 
sumed for stellar atmospheres. These can differ by up 
to 150-200 K, with higher temperatures resulting in es- 
timated Li abundances which are higher by ~ 0.08 dex 
per 100 K. Thus accounting for a difference of 0.5 dex 
between BBN and the observations, would require a se- 
rious offset of the stellar parameters. We note that there 
has been a recent analysis [l2| which does support higher 
temperatures, and brings the discrepancy between theory 
and observations to 2cr. 

Another potential source for systematic uncertainty 
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lies in the BBN calculation of the 7 Li abundance. The 
predictions for 7 Li carry the largest uncertainty of the 4 
light elements which stem from uncertainties in the nu- 
clear rates. The effect of changing the yields of certain 
BBN reactions was recently considered in Ref. @. In 
particular, they concentrated on the set of cross sections 
which affect 7 Li and are poorly determined both exper- 
imentally and theoretically. It was found for example, 
that an increase of the 7 Be(d,p)2 4 He reaction by a fac- 
tor of 100 would reduce the 7 Li abundance by a factor 
of about 3 in WMAP r] range. This reaction has since 
been remeasured and precludes this solution [13j ]. The 
possibility of systematic errors in the 3 He(a, 7) 7 Be reac- 
tion, which is the only important 7 Li production channel 
in BBN, was considered in detail in [14]]. However, the 
agreement between the standard solar model and solar 
neutrino data provides constraints on variations in this 
cross section. Using the standard solar model of Bah- 
call [IE], and recent solar neutrino data [HI], one can 
exclude systematic variations of the magnitude needed 
to resolve the BBN 7 Li problem at the > 95% CL QJ]. 
The "nuclear fix" to the 7 Li BBN problem is unlikely. 

On the other hand, various theoretical explanations 
involving physics beyond the standard model have been 
proposed |17|. One possible extension of the standard 
BBN scenario allows for inhomogeneous nucleosynthe- 
sis [l8[ but this seems to overproduce 7 Li. It has also been 
argued that particle decay after BBN could lower the 7 Li 
abundance and produce some 6 Li as well 19]. This has 
been investigated in the framework of the constrained 
minimal supersymmctric standard model if the lightest 
supcrsymmctric particle is assumed to be the gravitino 
[201 ]. Some models have been found which accomplish 
these goals [2l|. Another route is to assume that gravity 
is not described by general relativity but is attracted to- 
ward general relativity during the cosmic evolution [22| . 
BBN has been extensively studied in that scenario (see 
e.g. Ref. [23]). The effect of the modification of grav- 
ity is mainly to induce a time variation of the equivalent 
speed-up that can be tuned to happen during BBN but it 
can have other signatures both on cosmological and local 
scales [24| . 

In this article, we want to investigate the possible vari- 
ations of fundamental constants. It is well known that 
variations in the fundamental coupling constants such as 
the fine structure constant, a, can affect the light ele- 
ment abundance during BBN. Most analyses have con- 
centrated on the effect of such variations on the abun- 
dance of 4 He [i3 he m m m, m m m m m. 

Changes in a directly induce changes in the nucleon mass, 
Ani]v, which affects the neutron-to-proton ratio. 

Much of the recent excitement over the possibility of 
a time variation in the fine structure constant stems 
from a series of recent observations of quasar absorp- 
tion systems and a detailed fit of positions of the absorp- 
tion lines for several heavy elements using the "many- 
multiplet" method [34], l35fl . When this method is ap- 
plied to a set of Keck/Hires data, a statistically signif- 



icant trend for a variation in a was reported: Aa/a = 
(-0.54±0.12) x 10~ 5 over a redshift range 0.5 < z < 3.0. 
The minus sign indicates a smaller value of a in the past. 
In Ref. (36| , a set of high signal-to-noise systems yielded 
the result Aa/a = (-0.06 ± 0.06) x 10~ 5 . One should 
note that both results are sensitive to assumptions re- 
garding the isotopic abundances of elements used in the 
analysis which further complicates the interpretation of 
any positive signal from these analyses [37| . 

In addition to the possible variation in a, it is reason- 
able to search for other time-varying quantities such as 
the ratio of the proton-to-electron mass, /j, [H[. Indeed 
recent analyses [39j claim to observe a variation at the 
level 



Afi 



(2.4 ±0.6) x 10" 



(3) 



using Lyman bands of H2 spectra in two quasars. In the 
following, we investigate the effect on BBN of changes 
in fundamental couplings which could also account for 
a variation in fi. As we will see, the largest effect can 
be traced to a variation of the Higgs vacuum expecta- 
tion value leading to a variation in the binding energy 
of deuterium. We find that for a suitably large varia- 
tion (of order a few hundredths of a percent) in fi at the 
time of BBN, the 7 Li abundance can be decreased by the 
requisite factor without overly affecting the agreement 
between theory and observations for D and He. 

Although we do not directly tie our calculations to the 
observation with the result in Eq. , we do take as our 
starting point, the possibility that fi = m p /m e could have 
differed from its present value at the time of BBN. As a 
result we will be interested in variations of Yukawa cou- 
plings, h (we can assume, or not, that all Yukawas vary 
identically); the Higgs vacuum expectation value (vev), 
v, and A = h gcp . Some effects of the variations of v 
[30l . [3l| and A [32| on BBN have been considered in the 
past. Here we will be primarily interested in the inter- 
dependence of these variations (some of which are not 
model dependent) and their effects on quantities of di- 
rect importance to BBN, such as the binding energy of 
the deuteron, the nucleon mass difference and the neu- 
tron lifetime. 

In complete generality, the effect of varying constants 
on BBN predictions is difficult to model because of the 
intricate structure of QCD and its role in low energy 
nuclear reactions (see Refs. [IE HH). The abundances 
of light nuclei produced during BBN mainly depend on 
the value of a series of fundamental constants which in- 
clude, the gravitational constant G, the three gauge cou- 
plings and the Yukawa couplings of the electron and 
quarks. One needs to relate the nuclear parameters 
(cross-sections, binding energies and the masses of the 
light nuclei) to these fundamental constants. This ex- 
plains why most studies are restricted to a subset of con- 
stants, such as e.g. the gravitational constant 
the fine structure constant 0, HE [H, H3 or zTfscOiJ 

The approach we adopt here recognizes that many vari- 
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ations of parameters we deem fundamental are interre- 
lated in model dependent ways [13] • For example, if one 
assumes gauge coupling unification at some high energy 
(grand unified) scale, there is a direct relation between 
variations in the fine structure constant and A. In string 
theories, the variations of gauge couplings will be related 
to variations in Yukawa couplings, and in models where 
the weak scale is determined by dimensional transmu- 
tation [ill], there is a relation between variations in the 
Yukawa couplings and variations in the Higgs vev. Varia- 
tions in the latter will also trigger variations in A. While 
the exact relation between these variations is model de- 
pendent, the fact that they are interrelated is not. There- 
fore it is inconsistent for example to consider a variation 
in v without simultaneously varying A. We will make 
use of these dependencies to study variations in several 
(tractable) quantities which affect BBN. Coupled varia- 
tions of this type were used to strengthen existing bounds 
on the fine structure constant based on Oklo and mete- 
oritic data [44]. As noted above, we can not fully evolve 
the variations in all nuclear reactions, because their de- 
pendence is unknown. Here, we will be primarily inter- 
ested in the induced variations of the nucleon mass dif- 
ference, the neutron life-time, and the binding energy of 
deuterium. We recognize that this represents a limitation 
on our results. 

In Sec. UH we relate the variation of the BBN param- 
eters, mainly Q, t„ and Bd, to the variation of the fun- 
damental parameters such as the Yukawa couplings, h, 
the QCD energy scale, A, and the fine structure con- 
stant, a. Section HTT1 focuses on the relations that can be 
drawn between the variations of the fundamental param- 
eters, taking into account successively grand unification, 
dimensional transmutation and the possibility that the 
variation is driven by a dilaton. In order to deal with 
some of the theoretical uncertainties, we introduce two 
phenomenological parameters and we then make the con- 
nection with the variation of the proton to electron mass 
ratio at low redshift. Section IIVI focuses on the BBN 
computation and first describes the implementation of 
the variations in our BBN code. Assuming that the fine 
structure constant does not vary we show that deuterium 
and 4 He data set strong constraints on the variation of 
the Yukawa couplings [see Eq. (f3"9"| ] but that inside this 
bound there exists a range reconciling the 7 Li abundance 
with spectroscopic observations. We then allow the fine 
structure constant to vary and set a sharp constraint on 
its variation in the dilaton scenario [see Eqs. (|4"0")) and 

(EDI. 



FROM FUNDAMENTAL PARAMETERS TO 
BBN QUANTITIES 



II. 



As discussed above, we focus our attention on three 
physical quantities which have direct bearing on the re- 
sulting abundances from BBN, the nucleon mass differ- 
ence Q — m n — m„ = 1.29 MeV, the neutron lifetime r n , 



and the binding energy of deuterium Bd- 

The neutron-proton mass difference is expressed in 
terms of a, A, v, and the u and d quark Yukawa couplings 

as 



Q = m n — m p = a a A + (hd — h u )v . 



(4) 



where the electromagnetic contribution at present is 
a o;o Ao = —0.76, and therefore the weak contribution 
is (hdo — huo) vq = 2.05 [isl ]. The variation of Q will then 
scale as 



AQ 
Q 



= -0.6 



Aa 

a 



AA 



1.6 



A(h d -h u ) | At;' 
hd - h u v 



The neutron lifetime can be well approximated by 
_! 1 l + Sg\ 



(5) 



60 2tt 3 



■ Gl ml 



+ 151n(g+ y/q* - 1) 



(6) 



where q — Q/m e . Since Gp = l/\/2v 2 and m e = h e v we 
have for the relative variation of the neutron lifetime, 



Ar„ 



= -4., 



+3.8 



Av 
h 

v 

Aa 
a 



1.5 ■ 



AK 



10.4 



A(h d - h u ) 



AA\ 



(7) 



In addition to Q and r„, which have been well stud- 
ied in the context of BBN, we consider the variation of 
Bjj. This is one of the better known quantities in the 
nuclear domain: it is experimentally measured to a pre- 
cision better than 10~ 6 46] so that allowing a change 
of its value by a few % at BBN can only be reconciled 
with laboratory measurements if its value is varying with 
time. 

Recently, in a series of works [HI, HfJ, [4?J Flambaum 
and collaborators have considered the dependence of 
hadronic properties on quark masses and have set con- 
straints on the deuterium binding energy from BBN (STJ] 
following Refs. 0, HH HI HI . The importance of B D 
can be understood by the fact that the equilibrium abun- 
dance of deuterium and the reaction rate p(n, 7)D depend 
exponentially on Bjj and on the fact that the deuterium 
is in a shallow bound state. 

Here, we follow Refs. [2^, (4?J to compute the quark- 
mass dependence of the deuterium binding energy. Using 
a potential model, the dependence of -Bd on the nucleon 
mass and a and lu meson masses have been determined 



ABd 
B d 



50- 



6 



Ampf 
m N 



(8) 



for constant A. One can see that the coefficients of the 
quantities on the RHS of (jSJl do not add up to unity 
as it is required on dimensional grounds. Clearly there 
is a variation of a dimensional quantity that has not 
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been taken into account, which at low energy, can be 
expressed in terms of an a priori unknown combination 
of A and v. For definiteness we will write the missing 
term in terms of A only, keeping in mind that this some- 
what artificial method of fixing units is accounted for in 
the uncertainty in the relations between the variations of 
the fundamental parameters that we will discuss in the 
next section. Hence, when the nucleon and meson masses 
are kept constant, we write ABd/Bd = — 7AA/A. On 
the other hand, fixing A, when varying quark masses 
(the largest contribution comes from m s ), their result 
is ABd/Bd = —17 Am s /m s . 

The importance of the strange quark in the nucleon 
and meson masses can be traced to the 7r-nucleon £ term, 
which is given by 



ovjv = £ = - (m u 



m d )(B u +B d ). 



where 



Bn 



(p\qq\p) 



(9) 



(10) 



Defining y = 2B S /(B U + B d ), the combination £(1 — y) is 
the change in the nucleon mass due to the non-zero u, d 
quark masses, which is estimated on the basis of octet 
baryon mass differences to be Co = 36 ± 7 MeV [48[ . Fol- 
lowing Ref. H|, we have (B u - B s )/ \B d - B s ) = 1.49 and 
given a value of £, one can determine B q . In Ref. (47[, the 
value B s = 1.5 was adopted and corresponds to £ ~ 51 
MeV, which is a reasonable value. This corresponds to 



Ampf 
tun 



TUN 



0.19- 



For these values we find a similar value (though slightly 
larger than the one found in Ref. for the light quark 
(it and d) contributions which give 



Amjv 
tun 



0.052- 



This implies that 



Am, n 



AA [Ah 
0.76— + 0.24 — 
A V h 



Av 

v 



(11) 



The value of £ however has substantial uncertainties 
which were recently discussed in Ref. [50]| . A often used 
value is E = 45 MeV which was already somewhat larger 
than naive quark model estimates, and corresponded to 
y ~ 0.2. However, recent determinations of the 7r-nucleon 
E term have found higher values [5l| . E = 64 MeV. Still 
higher values can be ascertained for the observation of 
exotic baryons [52j|. For E = 45 (64) MeV, B s = 0.9 
(2.8) and Am N /m N = 0.12 (0.36) Am s /m s . The con- 
tribution from u and d quarks is 0.046 and 0.066 for E = 
45, 64 MeV, respectively. A similar calculation for the to 
meson leads to Am^/m^ = (0.09, 0.15, 0.29) Am s /m s 
for E = 45, 51, and 64 MeV respectively. For the a 



meson, three contributions were identified in [47]] , only 
one of which is related to E, yielding Ara a jra a — (0.44, 
0.54, 0.75) Am s /m s . Combining these sensitivities using 
Eq. ©, we would arrive at AB D /B D = (-16, -17, -19) 
Am s /m s (when the u and d contributions are neglected). 
Thus, despite the large uncertainties in the individual 
sensitivities, the dependence of Bd on the strange quark 
mass is relatively stable. Because of the cancellations 
in Eq. ([8]), the u and d quark contributions are indeed 
small: AB D /B D = (0.08,-0.009,-0.20) Am q /m q and 
can safely be neglected. 

Choosing the central value E = 51 MeV and since 
m s = h s v, we immediately have the relation between 
Bd, h, and v. Adding these two contributions and using 
ABd/Bd = —VjAm s /m s we have in general, 



AB 



D 



B D 



18- 



AA 



17 



Av Ah, 



(12) 



where once again we have repaired the mass dimension 
by adding the appropriate powers of A A/ A. 

Eqs. ([5]), (O, and (fT2]) form the initial basis for our 
computation. 



III. RELATIONS BETWEEN FUNDAMENTAL 
PARAMETERS 

A. General relations in a GUT 

We note that several relations among our fundamental 
parameters can be found. First, changes in either h or 
v trigger changes in A 53]. This is evident from the 
low energy expression for A when mass thresholds are 
included 



A = n 



m c m b m t 
M 3 



2/27 



exp 



2tt 



9a s (fi) 



(13) 



for /i > m t up to some unification scale in the standard 
model 1 



AA 



R 



Ac 



2 

27 



Ai 



Ah c Ahh 



h c 



hh 



Ah t 



The value of R is determined by the particular grand 
unified theory and particle content which control both 
the value of q:(Mgut) — & s (Mgut) an d the low energy 
relation between a and a s , leading to significant model 
dependence in R [54, 55] . Here we will assume a value of 
R = 36 corresponding to a set of minimal assumptions 
[27L [5^. However, in most BBN computations, we will 
neglect the variation in a and therefore the precise value 
of R chosen will not affect our conclusions. Nevertheless, 



In suporsymmetric models, additional thresholds related to 
squark and gluino masses would affect this relation l~)4| . 
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the relation between h, v and A is quite robust and has 
been neglected in most studies discussing the effect of 
varying v (or varying Gf) [30Ll3lj|. 

For the quantities we are interested in, we now have 



AB D 
B D 
AQ 
Q 



= -13 



Av 

v 



Ah 
~h 



18R- 



Ac 



a 



(15) 



l.: ) r^ + ^-0.6(l + J R)^, (16) 
v a I a 



Av 



Ah 
~h 



3.8(1 + R) 



A 



a 



a 



(17) 



where we have assumed that all Yukawa couplings vary 
identically, Ahi/hi = Ah/h. For clarity, we have writ- 
ten only rounded values of the coefficients, however, 
the numerical computation of the light element abun- 
dances uses the more precise values. We also recall that 
AGp/Gp = -2Av/v and Am e /m e = Ah/h + Av/v. 



B. Interrelations between fundamental parameters 

Secondly, in all models in which the weak scale is de- 
termined by dimensional transmutation, changes in the 
largest Yukawa coupling, ht, will trigger changes in v [43[. 
In such cases, the Higgs vev is derived from some unified 
mass scale (or the Planck scale) and can be written as 
(see Ref. |l|) 



Mp exp 



8n 2 c 
"ft 



(18) 



where c is a constant of order unity. Indeed, in su- 
persymmetric models with unification conditions such as 
the constrained minimal supersymmetric standard model 
there is in general a significant amount of sensitiv- 
ity to the Yukawa couplings and the top quark Yukawa 
in particular. This sensitivity can be quantified by a fine- 
tuning measure defined by [58j 



d In mw 
d In en 



(19) 



where mw is the mass of the W boson and can be sub- 
stituted with v. The are the input parameters of the 
supersymmetric model and include h t . In regions of the 
parameter space which provide a suitable dark matter 
candidate [59j| . the total sensitivity A = y/Y^i A'i typi- 
cally ranges from 100 - 400 for which the top quark con- 
tribution is in the range At = 80 — 250. In models where 
the neutralino is more massive, A may surpass 1000 and 
At may be as large as ~ 500. 

Clearly there is a considerable model dependence in 
the relation between Aw and Ah t - Here we assume a 
relatively central value obtained from Eq. (|18[) with c ~ 
ho — !• In this case we have 



Av 



Ah 



- 16n 2 c- „ 

v h A 



160- 



Ah 



(20) 



but in light of the model dependence, we will set 
Av _ Ah 

= o — — , 

v h 



(21) 



hence defining 5* = d In v/d In h ~ At and keeping in 
mind that 5 ~ 160. It follows that the variations of Bp, 
Q and r„ are expressed in the following way 



AB D 
B D 
AQ 



Ar„ 



= -17(5+1 
= 1.6(5+1) 



Ah 

Ah 
T ~ 



U8 
0.6 



4.85)4^ + 3. 
h 



AA 
"A - ' 

Aq A A 

a A 

Aa AA 
a A 



(22) 
(23) 
(24) 



where we have again assumed common variations in all 
of the Yukawa couplings. It also follows that AGf /Gf — 
-2SAh/h and Am e /m e = (1 + S)Ah/h. 
Now, using the relation (fl4|) we arrive at 



AB 



D 



B D 
AQ 

Q 
Ar„ 



-13(1 



-5) 
5) 

(8 + 45) 



1.5(1 



Ah 
~h 
Ah 
~h ~ 
Ah 
~ + 



18R 



Aa 



0.6(1 + R) 
3.8(1 + R) 



Aa 
a 
Aa 



(25) 
(26) 
(27) 



Finally we can take into account the possibility that 
the variation of the constants is induced by an evolv- 
ing dilaton [27|. In this scenario, it was shown that 
Ah/h = (l/2)Aa/a, therefore the expressions above can 
be simplified to 



^ = -[6.5{l + S)-im — 
Bd a 



AQ 

Q 

Ar„ 



(0.1 + 0.75 -0.6R) 
-[0.2 + 25-3.8J?] 



a 
Aa 



(28) 
(29) 
(30) 



though these relations will also be affected by model de- 
pendent threshold corrections. 



C. Sensitivity of Bd to the pion mass 

An independent calculation suggests a large depen- 
dence of the binding energy of the deuteron to the pion 
mass [13] parametrized in Ref. [3l|, for constant A, by 



AB 



D 



Am, 



Br 



(31) 



where r is a fitting parameter found to be between 6 
and 10. The mass of the pion is given by f^m 2 . = (m u + 
mci)(qq}, where f v oc A is a coupling and (qq) oc A 3 is the 
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quark condensate. Hence, the sensitivity of the binding 
energy to the fundamental parameters is 

^ = (> + i)^-i(v + x)' <*> 

which must be compared with Eq. (fT2|) . The coefficients 
are different, at most, by a factor of 4. Substituting 
relations (JTU and ([21]) into Eq. (J32J) we obtain 

^ = ( .2-0.4r)(l + S)^ + (l + 0.5r)E — , (33) 



which is to be compared with Eq. (f25j) . 

For the dilaton model considered, Ah/h = (l/2)Aa/a, 
and we find 

- [(0.1 - 0.2r)(l + 5) + (1 + 0.5r)i2] — . (34) 
Be, 1 a 

Again, the ratio of the coefficients in Eq. (|28|) to these is 
at the most of order 4. We therefore conclude that taking 
a scaled dependence with the pion mass (|3Tj) leads to 
the same constraints on the variation of the fundamental 
parameters (up to a factor of a few) as the nucleon/ meson 
mass dependence ©. 



D. Links to the variation of m p /m e 

Before we use the above relations in our BBN code, it 
is interesting to first compare these relations with the ob- 
served variation in fi. Using Eq. pip , and then Eqs. (TT4|) 
and (|2ip . the proton-to-electron mass ratio, /i = m p /m e 
varies according to 

^ = .8i?^- 0.6(5+1)^. (35) 
fx a h 

Using the current value on the observational variation of 
/i at redshift z ~ 3 [39], i.e. A/i/n w 3 x 10~ 5 we obtain, 
assuming a constant, 



Ah , / 161 

— ~ -3.2 x 10~ 7 

h \1 + S 

Interestingly we deduce from Eq. 
constant 



(36) 



5l) that when a is 



AB D 
B D 



22 



A^i 



6.6xl0~ 4 , (37) 



at z ~ 3, independent of the value of S. 

In the case where the variation is driven by a dilaton, 
we can link the observational variation in /i to a variation 
in a to get 



Aa 

a 



-1.5 x 10" 



-20.2 



0.8R - 0.3(5 + 1) 



(38) 



but higher than the stronger bound found in Ref. [36|, 
for the considered value (R,S) = (36,160). Note that 
this corresponds to ABd/Bd ~ 6 x 10~ 4 , by apply- 
ing Eq. I|28p. which is comparable to the value found in 
Eq. ([37]) where a was taken to be constant. 



which is compatible with the measurement of the time 
variation of the fine structure constant in Refs. 34, 35] 



IV. NUCLEOSYNTHESIS 
A. Implementation in a BBN code 

We incorporate the relations derived above in a BBN 
network. We use the reaction rates provided by the com- 
pilation and analysis of experimental data of Ref. (6l| 
covering ten of the twelve nuclear reactions involved in 
the BBN. The two remaining reactions of importance, 
n p and p(n, 7)D come from theory and are numer- 
ically evaluated, taking into account the variations dis- 
cussed above. 

The p(n,7)D reaction rate is calculated according to 
the Chen and Savage [62j derivation of the cross section in 
the framework of effective field theory. The weak reaction 
rates n <-> p are calculated by numerical integration of 
the electron, positron and neutrino Fermi distributions 
and phase space for the six weak interaction reactions. 
Zero temperature radiative corrections [63|, to the 
weak rates are also evaluated numerically. The neutrino 
versus photon temperature used in these rate calculations 
is a byproduct of the numerical integration of the Fried- 
mann equation where we follow the electron-positron an- 
nihilation exactly by the numerical integration of their 
Fermi distributions. 

In principle, one could have introduced a variation in 
Bd in the effective field theory cross section provided 
the scattering length in the 1 Sq channel follows this vari- 
ation. We found it simpler to use the Dmitriev et al. [I^ 
prescription for the reaction rate change that takes into 
account these effects. The binding energy of the deuteron 
is also directly involved in the calculation of the reverse 
rate. 

Variations in the binding energy of the dineutron also 
have little effect on the primordial abundances provided 
its absolute value remains smaller than the deuteron's 
binding energy [65]. Considering that, in this work the 
variations on the binding energy of the deuteron are only 
of a few percent, we do not expect any important role 
played by the binding energy of the dineutron in the cal- 
culations. 

By varying Bd one also changes the size of the 
deuteron which consequently modifies the D(D,n) 3 H c , 
D(D,p)T and D(p, 7) 3 H reaction rates [66]. We have, 
however, not included this effect in our calculation. 

The n <-> p weak rates depend on both Q and r n . 
Eq. ^ gives the decay rate of free neutrons; that is, the 
zero temperature limit of the weak n — -> p rate. This 
equation is used, in conjunction with the experimental 
value of the neutron lifetime (t„ = 885.7 s; see Ref. [671 ]) 
to fix the normalization of the finite temperature weak 
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rates. Hence, variations in r„ directly affect those rates 2 . 
More precisely, we use Eq. ([B]) to fix the present day nor- 
malization, then scale it according to the values of m e 
and Gf at BBN and use the BBN values of Q and m e to 
evaluate the integrals in the finite temperature weak rate 
calculations involving Q/m e , T/m e and T u /m e . Never- 
theless, the approximation of scaling the weak rates with 
T n is very good. 



B. Results under the hypothesis Aa/a — 

We proceed with the calculation of the light element 
abundances during nucleosy thesis using Eqs. (|25|) . ([26]) 
and (|27p first neglecting the contribution from Aa/a. 
From the discussions of the two previous sections, we 
conclude that variations in the Yukawa couplings lead 
directly to the changes in r„, Q, and Bjj. In addition, 
the nucleon mass also changes (Amjsr/rriN — 0.76AA/A+ 
.2A(Av / v + Ah / h) ~ 66 Ah/h) as does the electron mass 
(Am e /m e ~ l&lAh/h). 

In Fig. [TJ we show the time evolution of the light el- 
ements for the standard BBN model (solid curves) and 
when a variation Ah/h = 1.5 x 10 -5 (dashed curves) is 
considered, assuming the WMAP value for 77 and using 
the central value S = 160. We see that the largest ef- 
fect is indeed a decrease in the 7 Be abundance (which 
contributes to the final 7 Li abundance after decay) cor- 
related to an increase of n/p and a slight increase in the 
D/H abundance. 

In Fig. [21 we show the final abundances of D/H, 3 He, 
4 He, and 7 Li/H as a function of the variation in the 
Yukawa coupling h, for three assumed values of the pa- 
rameter S. The horizontal cross-hatched regions indicate 
the current observational spectroscopic determinations. 
For 4 He, we use Y = 0.232 - 0.258 70]. For D/H, wc 
use 2a range based the latest average of six quasar ab- 
sorption systems, D/H = (2.83 ± 0.52) x 10" 5 [zH - For 
7 Li/H, we show two ranges: the first given by 7 Li/H = 
1.23±g:|| x 10" 10 [I and the second given by 7 Li/H = 
(2.33 ± 0.64) x 10 -10 when higher surface temperatures 
are assumed fl2j and is represented with dashed lines. 
The dotted vertical line indicates the standard BBN re- 
sults (i.e. Ah/h = 0) for n = 6.12 x lO" 10 . 

We recall that there is significant model dependence in 
several of the assumed relations between the fundamen- 
tal parameters. For example, in Eq. (|21[) . we adopted 
c/h 2 = 1 (that is S ~ 160). However, the origin of 
the dependence between v and h depends on physics be- 
yond the standard model, and c/h 2 could be significantly 
larger or smaller than unity. 



2 While a more recent determination of neutron lifetime has been 
published [oil , its value has not been adopted by the RPP |67|1 . 
Consequences of this measurement on BBN were considered in 



Ah/h = and 1.5x10 




10 10 2 10 3 10 4 



Time (s) 

FIG. 1: The time evolution of the light element abundances 
for standard BBN (solid curves) and when variations Ah/h = 
1.5 x 10 -5 are included (dashed), assuming a remains con- 
stant. The WMAP value of r? = 6.12 x lCT 10 was assumed 
and we have used S = 160. 

As one can see from Fig. [2J each of the light elements, 
D, 4 He, and 7 Li show strong dependences on Ah/h. In 
fact, D/H provides us with the strongest constraint (un- 
der the hypothesis that a is constant) which for S = 160 
is, 

- 1.5 x 10~ 5 < y < 1.9 x 10~ 5 . (39) 

Using Eq. ([23]) . this bound translates to —4 x 10~ 2 < 
ABd/Bjj < 3.1 x 10~ 2 . Note that we have not used the 
7 Li abundance to set the lower bound on Ah/h. However, 
we also observe that for values of Ah/h > 1.8 x 10~ 5 
(Ah/h > 0.9 x 10~ 5 for the second range of observational 
7 Li), the 7 Li abundance is sufficiently small so as to come 
into agreement with the observational data. So long as 
we do not exceed the upper bound given in Eq. (|39|) . all 
of the light elements can be brought into agreement with 
data. Thus we must saturate the limit, but recall that 
this conclusion is based under the restrictive assumption 
that a is constant. On the other hand, it also means 
that our hypothesis can be falsified by decreasing the 
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S = 80, 160, 320, Aoc/oc=0 
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FIG. 2: Primordial abundances of 4 He, D, 3 He and 7 Li when 
allowing a variation of the Yukawa couplings. The horizontal 
cross-hatched regions depict the 2a spectroscopic data. We 
have assumed 3 values for S: S = 80 (dashed lines), S = 160 
(solid lines) and S = 320 (dot-dashed lines) for 77 = 6.12 x 

ht 10 . 



error bars of either 7 Li or D. 

In Fig. [31 we show the individual contributions of the 
varying BBN quantities to the light element abundances. 
When varying Q and m e individually, r« is held constant, 
i.e. Eqs. ([6]) and ([7]) have not been used. The effects of 
the variations of r„, Q, Bp, and m e can be seen ex- 
plicitly. The curves for m e are due to the effects of the 
electron mass on the expansion of the universe: m e ef- 
fectively enters in the r.h.s. of the Friedmann equation, 
affecting the timing and magnitude of the photon bath 
reheating following electron-positron annihilations. This 
effect is however very small as seen in Fig. [3] The effect 
of varying m e in the weak rates is accounted for in the 
overall variation of r„ . The electron mass does not affect 
the abundance of any of the isotopes, however, r„ and 
Q have a significant effect on 4 He leaving deuterium and 



7 Li almost unchanged. 

From the 4 He data, we deduce the bounds, —7.5 x 
1(T 2 < AB D /B D < 6.5 x 1CT 2 , -8.2 x 1CT 2 < At„/t„ < 
6 x 10" 2 and -4 x 1CT 2 < AQ/Q < 2.7 x 1(T 2 at 2a. A 
variation of the deuterium binding energy affects all the 
abundances, in particular, the deuterium data sets the 
tighter constraint -4 x 1CT 2 < AB D /B D < 3 x 1(T 2 . 
Interestingly, these bounds are equivalent to the ones ob- 
tained from the constraint (|39[) considering the interre- 
lations between the fundamental parameters. The 7 Li 
abundance is brought in concordance with spectroscopic 
observations provided its change falls within the interval 
-7.5 xl(T 2 < AB D /B D < -4xlCT 2 . We thus conclude 
that Bjj is the most important parameter connected to 
the discrepancy of the 7 Li abundance, and again, we see 
that there exists a window allowing for consistent 7 Li and 
deuterium abundances with data. 



m e , B D , Q and x n variations 



c 


0.26 


_o 








u 

C5 


0.25 












0.24 






i 






0.23 




0.22 




5 

a 



10 



m e 

— X 


D 






1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


3 He 

1 ! 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



10 





— — — - Qnp 


7 U^— 








mm\\\\\\\\\\\\\^ 


Mill 


1 1 1 1 1 1 1 1 1 1 1 1 1 


! 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



-0.1 -0.08-0.06-0.04-0.02 0.02 0.04 0.06 0.08 0.1 



Ax/x 



FIG. 3: Primordial abundances of the light nuclei as a func- 
tion of the relative variation of m e (dotted lines), r„ (dot- 
dashed lines), Q (dashed lines), and Bd (solid lines) with the 
same conventions as in Fig. [2] 

One may also consider the effect of the variation of 
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the nucleoli mass. The proton and neutron reduced mass 
enters as a factor (m" 1 + m" 1 )' in the p(n,j)D rate. 
For variations of the order we are considering, this effect 
is negligible. 



C. Allowing for Aa/a 7^ 

We now allow the fine structure constant to vary and 
we further assume that it is tied to the variation of the 
Yukawa couplings according to Ah/h — (1/2) Aa/a, us- 
ing Eqs. (j28|) - (|30|) . The results are shown in Fig. [4] where 
the abundances are depicted for three values of the pa- 
rameter R. Comparison of this figure with Fig. [2] shows 
the effect of including the variation in a. Not consider- 
ing 7 Li, the tighter bounds on Ah/h are again given by 
the deuterium abundance and are comparable in order of 
magnitude to the ones found in Eq. 



S = 240, R = 0, 36, 60, Aoc/oc=2Ah/h 



1.6 x icn 6 < 



Ah 
~h 



for R = 36 and 



, Ah 
- 3 x !Q-° < — 



< 2.1 x 10" 



< 4 x 10~ 5 , 



(40) 



(41) 



for S = 240 and R = 60. 

While these limits are far more stringent than the one 
found in Ref. [25|, it is consistent with those derived in 
Refs. [Ill [2?| where coupled variations were considered. 
Once again, for a variation near the upper end of the 
range (|40[) and (|4T|) . we can simultaneously fit all of the 
observed abundances. 

As noted above, a variation of a induces a multitude 
of changes in nuclear cross sections that have not been 
included here. We have checked, however, that a varia- 
tion of Aa/a w 4 x 10 -5 leads to variations in the reac- 
tion rates (numerically fit), mainly through the Coulomb 
barrier, of the most important a-dependent reactions in 
BBN [25] that never exceed one tenth of a percent in 
magnitude. 

Before concluding, we return once more to the ques- 
tion of model dependence. We have parametrized the 
uncertainty between Av and Ah with the quantity S and 
the uncertainty between AA and Aq through R. In full 
generality we ought to include one more unknown, say 
T, that parametrizes the relation between Aq and Ah, 
T = dhih/ dhia [5^|. In this work, however, we focused 
our investigation in the dilaton model where T — 1/2. It 
is now important to evaluate more precisely how sensitive 
our results are to the value these parameters may take. 
In Fig. [5] we illustrate the evolution of the primordial 
abundances of the light nuclei with S for a fixed value of 
the change in the Yukawa couplings assuming Aa/a = 0. 
We clearly see that, in this case, the theoretical 7 Li abun- 
dance is compatible with its observational measurement 
provided 200 < S < 370 (for the lower range of observa- 
tional Li abundances). 
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FIG. 4: Primordial abundances of 4 He, D, 3 He and 7 Li as a 
function of Ah/h — (1/2) Aa/a when allowing a variation of 
the fine structure constant for three values of the R parameter: 
R = (dashed lines), R — 36 (solid lines) and R — 60 (dot- 
dashed lines). 



We can also evaluate the impact of changing R in the 
dilaton model, when we allow a variation in a. To this 
end we show in Fig. [6] the evolution of the primordial 
abundances for two different values of Ah/h. We ob- 
serve that when Ah/h = 1.5 x 10~ 5 , we require R = 16. 
On the other hand, if we take Ah/h = 2.5 x 10~ 5 , the 
abundances are more sensitive to the value of R as the 
slope of the corresponding curves are steeper, but there 
is also a narrow window around R — 45 where all the 
light nuclei abundances are compatible with the full ob- 
servational data. 



V. SUMMARY 

In this article, we have considered the influence of 
a possible variation of the fundamental constants on 
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FIG. 5: Primordial abundances of the light nuclei as a func- 
tion of the parameter S assuming a change in the Yukawa 
couplings Ah/h = 1.5 x 10~ 5 and Aa/a = 0. 



FIG. 6: Primordial abundances of the light nuclei as a func- 
tion of the parameter R assuming a change in the Yukawa cou- 
plings: Ah/h = 1.5 X 1CT 5 (solid lines) and Ah/h = 2.5 x 1CT 5 
(dot-dashed lines), for S = 240. 



the abundances of the light elements synthesized during 
BBN. We have focused our attention on three fundamen- 
tal quantities central to BBN, namely Q, r„ and Bp, the 
variation of which was related to the one of the funda- 
mental constants. Specifying our theoretical framework 
we have reduced the fundamental constants to two inde- 
pendent ones, the Yukawa coupling and the fine structure 
constant. 

We have shown that these constants have a strong ef- 
fect on 4 He allowing us to set strong constraints on the 
variation of m e , Bd, Q and r n . Interestingly, the deu- 
terium and Li abundances are mainly sensitive to Bp 
and we have shown that there is a window in which 7 Li 
is compatible with spectroscopic data (see Fig. [3]). The 
existence of such a narrow window also implies that our 
mechanism can be falsified by an increase of the preci- 
sion of deuterium and/or 7 Li data. Our analysis also 
enables one to set sharper constraints on the variation of 



the fundamental constants. 

Assuming that the fine structure constant does not 
vary, we have shown that deuterium and 4 He data set 
strong constraints on the variation of the Yukawa cou- 
plings [see Eq. ([55)) ] but that inside this bound there 
exists a range reconciling the 7 Li abundance with spec- 
troscopic observations. We then allow the fine structure 
constant to vary and set a sharp constraint on its varia- 
tion in the dilaton scenario [see Eqs. p7J|) and l|4ip]. The 
theoretical limitations have also been discussed in detail. 
More specifically, we have parametrized the relations be- 
tween Aw and Ah and between AA and Aa with two 
free quantities, S and R, respectively. We found that the 
specific value of these quantities plays an important role 
alongside the change in the fundamental parameters in 
solving the 7 Li abundance problem. We conclude there- 
fore, that a better understanding of the values of these 
parameters from the theoretical standpoint can help us 
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to better constrain the variation of the fundamental pa- 
rameters at the time of BBN. 
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